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Abstract. Assembly line balancing problems consist in partitioning the 
work necessary to assemble a number of products among different sta- 
tions of an assembly line. We present a hybrid approach for solving such 
problems, which combines constraint programming and integer program- 
ming. 

1 Introduction 

Assembly lines are special flow-line production systems typical for the industrial 
production of high quantity standardized commodities. An assembly line consists 
of a number of work stations arranged along a conveyor belt. The work pieces are 
consecutively launched down the conveyor belt and are moved from one station 
to the next. At each station, one or several tasks necessary to manufacture the 
product are performed. The problem of partitioning the various tasks among 
the stations with respect to some objective function is called the assembly line 
balancing problem (ALBP) Jl6[ . 

Various classes of assembly line balancing problems have been studied in 
the literature. We will consider here so-called simple assembly line balancing 
problems SALBP of the following form: Let M = {1, . . . , m} be the set of stations 
and N = {1, . . . , n} be the set of tasks. We denote by tj the time required for task 
j and by Sj C M the set of stations able to fulfill task j. There is a precedence 
relation on the tasks that will be represented by a graph G = (N,E), where 
(iij 32) S E means that task ji is an immediate predecessor of task j 2 . By CTi 
we denote the total time available for executing the tasks that have been assigned 
to station i. The capacities CTi may vary from one station to another due to, 
e.g., different numbers of operators. 

The constraint satisfaction problem in assembly line balancing consists in 
assigning tasks to stations such that the total running time of the tasks assigned 
to some station does not exceed its capacity, and such that the precedence re- 
lations between the tasks are satisfied. The objective function is to minimize 
the number of stations necessary to fulfill all tasks. Including bin packing as a 
special case, the assembly line balancing problem is NP-hard. 
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The aim of this paper is to present a hybrid solver for assembly line balanc- 
ing problems, which combines constraint programming (CP) and integer pro- 
gramming (IP). The integration of integer programming and constraint pro- 
gramming has been an important research topic during the last years, see e.g. 
111^,000. The contribution of this paper is twofold: we develop a branch- 
and-cut solver for SALBP and show how it can cooperate with a CP solver in 
order to prune the search tree. 

The organization of the paper is as follows. We start in Sect. 2 with an inte- 
ger programming model of the simple assembly line balancing problem. Sect. 3 
describes the cutting planes that are used in the branch- and-cut solver on the 
IP side. Sect. 4 introduces a CP model of SALBP and Sect. 5 describes the co- 
operation between the IP and the CP solver. Finally, Sect. 6 contains a number 
of empirical results illustrating the benefits of the approach. 

2 Integer Programming Model 

Let Ti — {j : i € Sj} be the set of tasks which can be carried out by station i, 
and let A = Uj^i^j x {j}). We define the following decision variables 

{1, if task j is assigned to station i, 
0, otherwise. 

2.1 Constraints 

Feasible solutions x = [xij] of the assembly line balancing problem have to satisfy 
the following constraints: 

53 x ij = !. 3 G N > (1) 

53 tjXij < CT,, i G M, (2) 
jer, 

53 x ui ~ 53 Xi h - °' k — 1 , . . . , m : (ji,j 2 ) G E, (3) 

e {0,1}, {hj) G A. (4) 

SOS (Special Ordered Set) constraints ([!]) ensure that each task is assigned 
to exactly one workstation. Knapsack constraints (||) guarantee that the total 
running time of the tasks assigned to some station does not exceed its capacity. 
The constraints (||) correspond to the precedences given in the graph G. They 
express that if (ii, J2) G E and task j'2 is assigned to station k, then task ji must 
be assigned to one of the stations 1, . . . , k. 

2.2 Objective function 

Our objective is to minimize the number of stations necessary to perform all the 
tasks. Assuming w.l.o.g. that at most n stations are available, we define costs 



Cij for assigning task j to station i that satisfy the following condition: 

= Ci for all j G T t ; nc t < c i+ i, i = 1, . . . ,m — 1. (5) 

This ensures that lower numbered stations will be used first. Now the objective 
function can be written as follows: 

min ^2 E CiXi J " ( 6 ) 



3 Valid Inequalities and Cut Generation 

In order to solve the IP model presented above, we will use a branch- and- cut 
approach jl8| . In this section, we present the different classes of inequalities that 
will be used in our branch-and-cut algorithm. We start by defining the SALB 
polytope 

Psalb = convex G {0, 1} A : x satisfies d),©,©}). 

Here conv(S') denotes the convex hull of a set of points S C R n . We relaxed the 
SOS constraints (Q) to the inequalities 

x ij < 1. J6 N - (7) 

Since the polytope Psalb is contained in the multiple knapsack polytope (MK), 
defined as the convex hull of the set of points x G {0, 1} A satisfying (0) and (||), 
all inequalities valid for the MK polytope are also valid for the SALB polytope. 
Furthermore, the multiple knapsack problem is a special type of the generalized 
assignment problem (GAP); therefore, the inequalities valid for the GAP poly- 
tope are also valid for the SALB polytope. The known classes of inequalities for 
the GAP and MK polytopes are based on the notion of cover, which we briefly 
recall in the next section. 



3.1 Lifted Cover and (l,d)-Configuration Inequalities 

For a G JR" and b G M> , let 



P(a, b) = conv < x G {0, 1}™ : djXj < b 
{ J =1 

denote the knapsack polytope. A subset C C N is called a cover if X^ec a i ^ ^' 
a cover C is minimal if C \ {s} is not a cover for all s £ C. For a cover C, the 
cover inequality 

E x i < i^i - 1 ( 8 ) 



is valid for P(a, b); moreover, if C is a minimal cover, then (g) defines a facet of 

P c (a, b) = P{a, b)H{xe R" : x 3 = 0, j £ N \ C}. 

A pair (H, z) is called a (1, d)- configuration |t2| if H C N, z € N \ H, and 
2 < d < \H\ are such that 

— H' U {z} is a minimal cover for every H' C H with = d. 
If (H,z) is a (1, d)-configuration, the inequality 

^ Xj- + (|JT| - d + l)x z < \H\ 
■jeH 

defines a facet of P Hu ^(a, b). 

Lifting of inequalities is a key issue in branch-and-cut. It allows one to 
strengthen an inequality by calculating non-zero coefficients for variables that 
initially are not present. For a formal definition and general results about lifting 
see Q. The lifted cover and (l,d)-configuration inequalities are automatically 
generated by the solver we have used. For the precedence-constrained knapsack 
polytope, a straightforward generalization of minimal covers, the so-called min- 
imal induced covers, has been investigated, see e.g. ||[l3|,[E7|. By analogy, we 
introduce the induced cover inequalities for the SALB polytope. 

Two tasks ji , j2 S N are called incomparable if both (ji , ji ) ^ E and ( j'2 2" 
E. A set W C N is called incomparable if the elements in W are pairwise 
incomparable. 

To each station k £ M corresponds a knapsack given by the inequality 
12j£T tjXkj < CTfe. We say that C C Tfc is a minimal induced cover (MIC) 

i if k 

— C is incomparable, 
~~ EjeC< fy^kj > CT k , 

~ Ejec<\{ s } ^ CT k for a11 s € C. 

Here C- d = {j <E N : j < ji for some j\ £ C}. 

Let C C TV be a MIC for knapsack k. Then the inequality 

jec j£C^\CieSj -.i<k 
is valid for the SALB polytope. 

3.2 Cycle Inequalities 

The class of cycle inequalities was introduced in § for the GAP polytope. Since 
we are not aware of a separation algorithm for this class, we describe here a 
separation heuristic for the subclass of cycle inequalities with cycles of length 4. 



Given two tasks u and v; assume that t u < t v . Let Ck and Ci be covers 
of the knapsacks k and I resp. such that: a) u,v S Cj,, b) it <E Cj, w ^ C/, c) 
C fc flC, = {u}. Then 

E *ii<|Cfcl + |tfl|-3 (9) 

ieCfe ieciu{ti} 

is a valid inequality for the SALB polytope || . 



4- Cycle Heuristic 

— Choose two knapsacks k, I 6 M, and two items it, u S 2fe n T; such that, 
for the point x = [xij] to be separated, all four values Xk u , %iu, %kv, %iv are 
greater than zero. Assume that t u <t v . 

— Compute a minimal cover C' k of the knapsack X)jer fc \{u v} tj x lj — ^Ti — 
t u - t v by solving 

min V (1 - x kj ) : V t, > CT k - t u -t v + l. (10) 

- * u 'jec jec 

set c fe = q; U {u, v}. 

— Compute a minimal cover C[ of the knapsack X^eT,\Cfc ^i^'i — ^^fe — ^« by 
solving 

nPfrn " Xy) : E*^ CTfc -tu + l. (11) 

Set Ci = C[ U {u}. 

— If both covers Ck and C\ exist, lift inequality (||) using the procedure de- 
scribed in section |3.4{ otherwise, return failure. 



3.3 Extended Cover and Heterogeneous Two-Cover Inequalities 

Let C be a cover with respect to some knapsack k G M, and £) C AT \ C be such 
that D U {i} is a cover of knapsack I € M\ {k}, for all i g C Then the extended 
cover inequality 

E(^ + + E ^ ^ i c ' + \°\ - 1 

J'EC iGD 

is valid for the multiple knapsack polytope ||. 

Consider two knapsacks k, I £ M, k I. Assume that C C Tk is a cover with 
respect to knapsack k, and let D be a subset of TJ \ C such that for all D' C D 
and C" C C with |C"j = \D'\, the set (C \ C) U D' is a cover for knapsack Z. 
Then the heterogeneous two- cover inequality 

5>jw+ E (|C|-i)^<|c|(|C|-i) 

is valid for the multiple knapsack polytope || . For separation of extended cover 
and heterogeneous two-cover inequalities, we use the heuristics described in B. 



3.4 A General Lifting Procedure 

Let A C A and _^ 

E a v x H<0 (12) 
(i.i)eA 

be a valid inequality for Psalb- For (i , j ) e A\A, the inequality 

if valid for Psalb if the coefficient cti j is computed by the following procedure: 
Lifting Heuristic 

— Let A* — A. For (i, j) G A, if j ~ jo, or io < i and (j, jo) G P, or io > i and 
(io>i) S P, exclude from A*. 

— For each i G M, 

• generate a permutation 7T 1 = (ir\, . . . ,7r^.) of the elements of Ti, = |Tj|; 

• for s = 2, . . . , n,-, compute 



^ max j]T> : y G {0, 1} S , J^ife < CTi j ■ ( 14 ) 

— Solve the optimization problem 

7 = max ^ onjZij 

E z « ^ i' j e ^> 

E 2 w* < < s - s = 2, . . . , m; i G M, 

l<j<s, 
(l,7T*)GA* 

2f tf G {0,1}, (t,j) G A*. 
and set o, j = j3 — 7. 

Note that both optimization problems (E_4|) and ( |l5| ) can be solved in polynomial 
time. Problem (|l4|) is a special type of knapsack problem and can be solved effi- 
ciently by dynamic programming. Problem (|l^) is a special case of the weighted 
matroid intersection problem . In fact, it can also be reduced to the minimum 
cost maximum flow problem, see e.g. H. 



4 Constraint Programming Model 



Our CP formulation of the SALB problem is based on the cumulative constraint 
[El. The basic version of cumulative can be defined as follows. 



There arc n tasks; task j is characterized by three parameters, which can be 
either domain variables or values: the starting time start j, the duration durj, 
and the amount resj of some resource consumed by the task. We are also given 
the completion time e for all the tasks, and the upper bound v on the resource 
consumption; e and v again are domain variables or values. The global constraint 

cumulative([[storii, dur\, res{\, . . . , [start n , dur n , res n ]], v, e) 

is satisfied if the following conditions hold: 

resj <v, t = 1, . . . , e, 

l<j<n : 
start j <t<startj +durj 

max [start j + durA < e . 

l<j<n 

In addition to the n given tasks, we introduce m artificial tasks numbered 
n+ 1, ...,n + m d = n; for CT max = max ieM CT h let S n+i = {i}, t n+l = 
CT max - CTi, i s M. Each task j G {1, . . . , n + m} is associated with a triple 
(start j,durj, res j) of domain variables, where 

— start j = i if task j is assigned to station i, 

— rfurj = 1, 

— res.j =tj. 

Now the CP model can be stated as follows: 

{ start j £ Sj, j = 1 , . . . , n, 
durj £ {1}, j = 1, . . . , n, 
resj £ {t,-}, j = 1, . . . ,n. 

!cumulative([[siarti, dttn, resi], . . . , [siarijj, durn, resn]], 
CT max ,m+l); 
start j! < start j 2 , (ji,j2) £ 

5 Combining IP and CP 
5.1 Reducing the Problem Size 

The size of the IP and CP model described before can be reduced as follows. 

In a first step, we build the CP model and do propagation, but without 
labeling. Let Sj denotes the domain of start j after propagation, and let Ti = 
{j : i £ Sj}. The IP model is obtained by imposing the constraints (0)-(|J) for 
these reduced sets Sj and Tj. 

In a second step, we do propagation by combining both solvers. Iteratively, 
for each task j £ N, we perform the following operations: 

— Minimize J2ieS- i' x ij subject to the constraints (|l|)-(|^) using cut generation, 
but without branching. Let 71 denote the objective value returned by the 
solver. For i £ Sj and i < [71 J , set Xij = and remove i from the domain 
of start j. 



— Maximize X^ies * ' Xi i un der the constraints using cut generation, 
but without branching. Let 72 denote the objective value returned by the 
solver. For i £ Sj and i > [72] , set Xij = and remove i from the domain 
of start j. 

— Initiate propagation for the CP problem. For q £ N, if i is not in the domain 
of start q , set Xi q = 0. 

5.2 Propagation and "Rounding Off' LP solutions 

After reducing the problem size, we can use both solvers, IP or CP, to continue. 
If we decide to use the CP solver, there is no further interaction with the IP 
solver. 

If we use the IP solver, we propose the following cooperation between IP and 
CP. Assume that Sj, for j £ N, is the domain of start j after preprocessing. For 
these sets Sj, consider the IP problem 

Propagation for subproblems. The IP procedure starts processing a node of 
the branch-and-cut tree by calling the CP solver in the following way. For the 
LP subproblem at this node, let d~j,dfj £ {0, 1} denote the current lower and 
upper bounds of the variables x^. For j £ N, set Sj = {i : d^- = 1}. For these 
sets Sj, build the CP problem and do propagation. At the end, let dfj = for 
those i £ Sj that are not in the domain of start j . 

Rounding of LP solutions. Let x — be a solution of some LP subprob- 
lem. For j £ N , define Sj = {i : Xij > 0}, build the CP model, and run labeling 
for a limited amount of time, trying to find a feasible solution. 

6 Empirical Results 

This paper describes ongoing work. Therefore, we can give only preliminary 
empirical results. These are based on a selection of benchmarks from the data 
sets in [ fT6| . For these examples, the cycle time CT is the same for all stations, i.e., 
CT = CTi,i £ M. The computational experiments were done on a Pentium III 
600 MHz using the CHIP C Library j| and the branch-and-cut code developed 
by the second author. 

We co mpare two versions of the problem reduction heuristics described in 
Sect, fy] The first version generates cuts only when computing an initial fea- 
sible LP solution. All other LP problems are solved without generating cuts 
(except standard lifted cover and (1, d)-configuration cuts, which are automati- 
cally produced by the IP solver). The second version generates all the different 
cuts for all LP problems. 

We measure the quality of a problem reduction strategy by the sum of the 
cardinalities of the sets Sj upon termination of the procedure. These values are 
given in the last three columns of Table [j] for CP alone, CP and LP with standard 
cuts, and finally for CP and LP with all cuts. 



Instance 
name 


# tasks 
n 


Cycle 
time CT 


# stations 
m 


Initial 
size 


After 
CP 


After LP with 


standard cuts 


all cuts 


Sawyer30 


30 


47 


7 


210 


112 


109 


33 


Sawyer30 


30 


28 


12 


360 


187 


187 


179 


Gunther 


35 


54 


9 


315 


129 


105 


105 


Gunther 


35 


44 


12 


420 


182 


182 


176 


Lutz3 


89 


118 


14 


1246 


285 


262 


195 


Lutz3 


89 


74 


23 


2047 


386 


386 


318 


Warnecke 


58 


155 


10 


580 


248 


241 


212 


Warnecke 


58 


73 


22 


1276 


575 


575 


575 


Tonge70 


70 


251 


14 


980 


394 


394 


394 



Table 1. Reducing problem size 



Next we give in Table |2j the total running times for the hybrid solver described 
in Sect. |5.2| including problem reduction by CP propagation. 



Instance 


CP 


Cut 


Total 


name 


prop. 


gener. 


sol. 




time 


time 


time 


Sawyer30 (47-7) 


1.261 


2.535 


5.809 


Sawyer30 (28-12) 


0.260 


0.901 


2.824 


Gunther (54-9) 


0.031 


0.230 


0.401 


Gunther (44-12) 


0.191 


0.941 


1.712 


Lutz3 (118-14) 


0.821 


14.151 


16.183 


Lutz3 (74-23) 


1.895 


42.019 


47.849 


Warnecke (155-10) 


5.597 


17.367 


38.345 


Warnecke (73-22) 


2.765 


29.007 


1:36.499 


Tonge70 (251-14) 


1:27.120 


10:22.7 


13:49.330 



Table 2. Running time in min:sec 



Note that neither the CP nor the IP solver alone are able to solve, for example, 
problem Warnecke (73-22) in less than one hour of running time. 

7 Conclusion and Further Research 

Current computational experience does not allow us to tell whether CP and IP, 
alone or in cooperation, are able to solve to optimality the unsolved instances 
of the simple assembly line balancing problem in the data sets from |Q. Fur- 
ther research and implementation work is necessary to answer this question, in 
particular when the precedence relations are not tight. 



While this is an interesting theoretical question, the main impact of this 
work is practical. The hybrid IP/CP solver described in this paper may become 
a platform for modeling and solving assembly line problems with all kind of side 
constraints, which are typical for real-world industrial applications. 
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